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Abstract: In this talk we present our detail study ( theory and numbers ) on the shad- 
owing corrections to the gluon structure functions for nuclei. Starting from rather contraversial 
information on the nucleon structure function which is originated by the recent HERA data, 
we develop the Glauber approach for the gluon density in a nucleus based on Mueller formula 
[H and estimate the value of the shadowing corrections in this case. Than we calculate the 
first corrections to the Glauber approach and show that these corrections are big. Based on 
this practical observation we suggest the new evolution equation which takes into account the 
shadowing corrections and solve it. We hope to convince you that the new evolution equation 
gives a good theoretical tool to treat the shadowing corrections for the gluons density in a 
nucleus and, therefore, it is able to provide the theoretically reliable initial conditions for the 
time evolution of the nucleus - nucleus cascade. The initial conditions should be fixed both 
theoretically and phenomenologically before to attack such more complicated problems as the 
mixture of hard and soft processes in nucleus-nucleus interactions at high energy or the theo- 
retically reliable approach to hadron or /and parton cascades for high energy nucleus- nucleus 
interaction. 



* E-mail: ayala@if.ufrgs.br 
** E-mail:gay@if.ufrgs.br 

t E-mail: levin@hep.anl.gov;leving(S)ccsg. tau.ac.il 



1 Introduction. 



The main goal of this talk is to share with you our experience and results that we got during 
the last two years reconsidering the whole issue of the shadowing corrections ( SC ) to the gluon 
density in nuclei [0]. The title which reflects the key problems that we are going to discuss is: 
"All ( theory and numbers ) about the SC to gluon density in nuclei" 

It is well known that the gluon density is the most important physical observable that 
governs the physics at high energy (low Bjorken x) in deep inelastic processes 0. Dealing with 
nucleus we have to take into account the shadowing correction (SC) due to rescattering of the 
gluon inside the nucleus, which is the main point of interest in this paper. We show that SC can 
be treated theoretically in the framework of perturbative QCD (pQCD) and can be calculated 
using the information on the behavior of the gluon structure function for the nucleon. 

The outline of the talk looks as follows. We start with our motivation answering the question 
why we got interested in the SC for nucleus gluon density. In section 3 we will discuss the 
theory and numerics of the Glauber ( Mueller ) approach emphasizing it's theory status and 
the estimates for the SC that came out of it. After short discussion in section 4 the first 
corrections to the Glauber approach we will present what we consider as a right way of doing, 
namely, the new evolution equation that sums all SC ( section 5 ). In section 6 we are going to 
discuss our next steps that we plan to do in a nearest future, while in section 7 we will give our 
answer to the hot question: and what ?, trying to collect all problems of RHIC physics that we 
will be able to answer using our approach. 



2 Motivation. 

Let us start with a brief summary of the HERA results for the nucleon structure functions 
( parton densities in a nucleon). The experiment |Q shows that the deep inelastic structure 
function F2{x, Q^) increases in the region of small x ( at high energies): 

F2(a;,Q2) oc for 10"^ > x > 10"^ 

at large and small ( ~ 1 — 2 GeV^ ) values of the photon virtualities Q^. 

At first sight we can conclude from the analysis based on the DGLAP evolution equations 
undertaken through all the world [H that : 

1. The DGLAP evolution equations work quite well and no other ingredients are needed to 
describe all the HERA data. 
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2. The parton cascade is rather deluted system of partons with small parton - parton 
interaction which can be neglected in a first approximation. In other words we do not need any 
SC to describe the experimental data. 

3. The phenomenological input, namely, the quark and gluon distribution at initial virtuality 
= Ql can be chosen at sufficiently low values of using the backward evolution of the 

experimental data in the region of ~ 4 — ^GeV^. Even more, the craziest parameterization 

that we have seen in our life - the GRV one does it's job perfectly well, starting with 
Ql = 0.3Ge\/2 ?, 

What we have discussed is moreless common opinion of all experts in DIS and one can find 
it in many plenary and review talks during the last two years. 

However we would like to draw your attention to several facts which do not fit to this 
common scheme: 

1. The best parameterization of the HERA data is not the solution of the DGLAP equations 
but a simple formula H]: 



F2(x, Q^) = a + m log ^ log — 



with a = 0.078 ; m = 0.364 ; Xo = 0.074 ; Ql =0.5 GeV^. It is clear that this simple formula 
cannot be a solution of the DGLAP evolution equations. To make obvious this remark it is 
enough to recalculate the gluon structure function from the above expression as it has been 
done in Ref.|P]. Indeed, xG^x.Q"^) turns out to be equal to 

xGix.Q'') = 3 log — 

X 

without any - dependence within the direct contradiction with the DGLAP evolution. 

2. Using the HERA data we can evaluate the parameter which characterizes the value of 
the SC, namely [§] 

^ = xG{x,Q'^) , (1) 

where xG{x, Q"^) is the gluon structure function and is the radius of area populated by 
gluons in a nucleon. The physical meaning of k becomes clear if we rewrite it in the form 



n = xG{x, Q 



2, a{GG) 



where a{GG) is the cross section of two gluon interaction in our parton cascade calculated by 
Mueller and Qiu |10|, namely, a{GG) = The physical meaning of this formula is the 

probability of the gluon - gluon interaction inside the parton cascade. It looks very natural 
if we compare eq. (|I]) with the small parameter for proton - nucleus interaction. Indeed, the 
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Figure 1: The J/\E' production without a) and with b) dissociation of the proton. 

parameter which governs the value of the Glauber corrections for proton - nucleus interaction 
KpA = A^^^, where A is the number of the constituents (nucleons) , a{pp) is the cross section 

of the interaction of our constituents and vri?^ is the area populated by nucleons. The question 
arises what is the value of R in eq. (|l|)? Using the new HERA data on photoproduction of 
J/\l/ meson [jlT| we are able to estimate the value of in the definition of k, (see eq. (|1|) ). 
To illustrate the point we picture in Fig.|I| the process of J/\l/ photoproduction in the additive 
quark model (AQM ). We see that we have two processes with different slopes {B ) in t { 
or in b]_ ): the J/\l/ production without ( Fig.|l]a ) {Bgi = bGeV'"^) and with ( Fig.|I]b ) ( 
Bin = 1-Q6 GeV~^ ) dissociation of the proton. The AQM gives us the simplest estimates for 
the resulting slope { R^ ) in eq. (|1|) if we neglect any slope from the Pomeron - J/\l/ vertex in 
FigJ^, namely 

Fig.^ shows the contour plot for k using the GRV parameterization P| for the gluon structure 
function and the value of R^ = 5GeV~^. One can see that k reaches k = 1 at HERA kinematic 
region, meaning shadowing corrections take place. 

3. The situation looks even more contraversparameterizationial if we plot the average value 
of the anomalous dimension < 7 > = dln{xG{x,Q'^))/dlnQ'^ in the GRV parameterization. 
1^ Fig. 3 shows two remarkable lines: < 7 >= 1, where the deep inelastic cross section reaches 
the value compatible with the geometrical size of the proton, and < 7 > = 1/2, which is the 
characteristic line in whose vicinity both the BFKL Pomeron ( see Ref.||12|| ) and the GLR 
equation should take over the DGLAP evolution equations. We will discuss later what are 

* We will discuss below the definition of the anomalous dimension and why this ratio is the average anomalous 
dimension. 
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Figure 2: Contour plot for k, for = 5Gef 

the BFKL and the GLR equations, what we need to know right now, is only the fact that both 
equations give the signal of the new physics. The HERA data passed over the second line and 
even for sufficiently small values of they crossed the first one without any indication of a 
strange behaviour near these lines. 

Concluding this brief summary of the HERA data and physics behind them we would like 
to repeat that to our taste the situation at HERA looks very controversial and the statement 
that the DGLAP evolution works is first but not the last outcome of the HERA data. On the 
other hand we have to develop the new approach to the SC, more general than the GLR one, 
which will allow us to give reliable estimates for the SC in the kinematic region to the left of 
the line < 7 > = 1/2. This is why we decided to reconsider everything that has been known 
about the SC, trying to forget everything that we knew about them, and to start our analisys 
of the SC from the very beginning. We also decide to choose the gluon density in a nucleus as 
a laboratory or training ground for the new approach to the problem of the SC. 

We have three reasons for such a choice: (i) the nucleus DIS is easier to handle theoretically, 
as we will show in the main body of our talk;(ii) the previous analysis of the SC shows that this 
is mainly density effect in the parton cascade ( see review ||13| for example ) and we anticipate 
larger gluon density for DIS with a nucleus; (iii) the RHIC is coming and the gluon density in 
nuclei will provide the initial condition for any phenomenological cascades for nucleus - nucleus 
interaction at high energies. 



3 The Glauber approach in QCD 

3. 1 The Mueller formula. 
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Figure 3: Contours for < j > = 1 and < 7 = 1/2. 




Figure 4: The structure of the parton cascade in the Glauber ( Mueller) formula. A denotes the 
nucleus, N - the nucleon, G*{Q^) - the virtual gluon and (7i^{r^) is the nucleon cross section. 
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The idea how to write the Glauber formula in QCD was originally formulated in two papers 



Ref . and Ref . . However, the key paper for our problem is the second paper of A. Mueller 
who considered the Glauber approach for the gluon structure function. Nevertheless, it is easier 
to explain the main idea considering the penetration of quark - antiquark pair, produced by 
the virtual photon, through the target. While the boson projectile is traversing the target, 
the distance r± between the quark and anti-quark can vary by amount Ar± oc Ra^ where E 
denotes the energy of the pair in the target rest frame and Ra is the size of the target (see 
Fig.^). The quark transverse momentum is kt oc l/r±. Therefore 

A r± oc < , (3) 

and is valid if 

r^s > 2mi?, (4) 
where s = 2mE. In terms of Bjorken x, the above condition looks as follows 

X < (5) 

Therefore the transverse distance between quark and antiquark is a good degree of freedom 
13 HOI- been shown by A.Mueller, not only quark - antiquark pairs can be considered 



in a such way. The propagation of a gluon through the target can be treated in a similar way 
as the interaction of gluon - gluon pair with definite transverse separation rt with the target. It 
is easy to understand if we remember that virtual colorless graviton or Higgs boson is a probe 
of the gluon density. 

The total cross section of the absorption of gluon(G'*) with virtuality and Bjorken x can 
be written in the form: 



a 



= fo'^' /l^vl/f (Q2^r„x,z)a^(x,r2) [vl/f (Q^ r„ x, z)]* , (6) 



tot^ 



where z is the fraction of energy which is carried by the gluon, \E'j^* is the wave function of the 
transverse polarized gluon and aA{x,r^) is the cross section of the interaction of the GG- pair 
with transverse separation with the nucleus. This cross section can be written in the form: 

a,(.,ri) = 2 I d\l,na(.,r,M , (7) 

where a is the elastic amplitude for which we have the s-channel unitarity constraint: 

2 Ima{x,r±,bt) = \a{x,r±,bt)\'^ + Gin{x,r^,bt) , (8) 

where Gm is the contribution of all the inelastic processes. Let us recall that two terms in 
eq. (M) have different physical meaning: the left hand side and the first term in the right hand 
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side describe the interference between the incoming plane wave and outgoing spherical wave 
which amplitude is the elastic scattering amplitude (a). These two terms cannot be calculated 
using a classical approach or simple Monte Carlo - like model. The Quantum Mechanics 
of the interaction is mostly absorbed in these two terms while the last term has a simple 
probabilistic meaning, namely, the probability of any inelastic interactions, and can be treated 
almost classically and, for certain, in the probabilistic way, for example in Monte Carlo-like 
models. The unitarity establishes the correlation between two unknowns a and Gin and has 
the general solution: 

a{x,r^A) = i{l - e-5f^(^'.'-i'''*) } ; (9) 

One can see that Vt has a simple physical meaning, namely is the probability that GG-pair 
has no inelastic interaction during the passage through the target. The opacity Vt is an arbitrary 
real function, which can be specified only in more detail theory or model approach than the 
unitarity constraint. One of such specific model is Glauber approach or Eikonal model. 

However, before we will discuss this model let us make one important remark on the strategy 
of the approach to the SC. We are trying to built a model or theory for the total cross section 
( or for the gluon structure function ) not because the SC should be the strongest one in 
this particular observable, but because if we will be able to calculate opacity VL we will have 
the theory or model for all inelastic processes. Indeed, using AGK cutting rules |jl5| we can 
calculated any inelastic process, if we know f2, in accordance with the s-channel unitarity. It 
is worthwhile mentioning that the inverse procedure does not work. If we know the SC in all 
details for a particular inelastic process, say for the inclusive production, we cannot reconstruct 
all other process and the total cross section in particular. 

Now, let us built the Glauber approach. First, let us assume that f2 is small (i7 ^ 1 
) and it's ht dependence can be factorised as i7 = VL{x,r^_)S{bt) with the normalization: 
/ d%tS{ht) = 1. Expanding eq. (^) and substituting it in eq. (|^, one can obtain: 

(^A{x,r±) = n{x,r±) (10) 

At small Q the cross section of the deep inelastic process with a nucleus is proportional to the 
number of nucleons in a nucleus {A) , namely. 

To calculate Q we need to substitute everything in eq. (^) and use the formula for cr^(G*) = 
^xGa{x,Q'^) as well as the expression for the wave function of the GG- pair in the virtual 
gluon probe. Such calculations has been done in Ref.0 and we recapture here the result ( see 
for example Ref. ^ for more details ): 

= AaNix,r±) = r']_xG{x, —) . (11) 
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The Glauber (eikonal ) approach is the assumption that = S{bt) with of eq. (|TT 



not only in the kinematic region where is small but everywhere. From the point of view 
of the structure of the final state this assumption means that the rich typical inelastic event 
was modeled as a sum of the diffraction dissociation of GG - pair plus uniform in rapidity 
distribution of produced gluons. For example, we neglected in the Glauber approach all rich 
structure of the large rapidity gap events including the diffractive dissociation in the region of 
large mass. 

Substituting everything in eq. ([Tl|) and eq. ® and using the wave function calculated by 
Mueller in Ref. we obtain the Glauber (Mueller) formula for the gluon structure function: 

xGAx,Q^) = ^j'^ r^^2{l-e-^^^^(^>?)-(^?)} (12) 

TT Jx X J TlTi Jo TT ^ ^ 



It is easy to see that the first term in the expansion of eq. (|T^) with respect to a gives the 
DGLAP equation in the region of small x. 

To calculate the profile function S{bt) we make the usual assumption that < bf >n in the 
interaction of GG - pair with the nucleon is much smaller than the nucleus radius ( < 6^ > <^ 
Ra- Therefore, S{bt) can be expressed through the nucleon wave function in a nucleus, namely 

A 

SA{qzA) = / dzi e'^^''<ilAizi,bt;r2,...ri,rAW{zubt]r2,...ri,rA) Y[ d^ri , (13) 
where the wave function is normalized as 

A 

/ <ilAizi,bt;r2,...ri,rAWizi,bt;r2,...ri,rA) H d^n = A. (14) 

i=i 

Assuming that there is no correlation between nucleons in a nucleus and the simple Gaussian 
form of a single nucleon wave wave function we derive the Gaussian parameterization for S{bt), 
namely 



j2 fj'2 

SAiq.X) = -^e~^A-^'\ (15) 

TT 



where the mean radius Ra is equal to 



2 

p2 _ ^ p2 
^A — g ^ws 



and Rws is the size of the nucleus in the Wood-Saxon parameterization. We choose Rws = 
To Aa with ro = 1.3 fm in all our calculation. We are doing all calculation in the rest frame of 
the nucleus where we can neglect the change of energy for the recoil nucleon in the nonrelativistic 
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theory for the nucleus. Indeed, its energy is Ep' = + ^ and ^ <^ g^. At high energy ( 
small X ) we can neglect also the g^-dependance ( see Ref.[^ for details). 

Using Gaussian parameterization for S{bt) ( see eq. (|T5|) ) we can take the integral over bt 
and obtain the answer {Nc = Nf = 3) 

xGa(x,Q^) = ^/'^ /^^{C' + M'^G(x',r,2))i?i(«:c(x',r,2))} (16) 

where C is the Euler constant and Ei is the exponential integral (see Ref.[^] Eq. 5.7.11) and 

The eq. (|T^ is the master equation of this section and it gives a way to estimate the value of 
the SC. We would like to stress that we have only adjusted the approach of Ref. to the 
rescattering in a nucleus. It means that we did nothing except that we share the responsibility 
with A. Mueller for eq. (||). 

One can see that the Mueller formula of eq. (|T2D depends only on k. If n is small ( /t ^ 1 
), we can expand eq. (|T6D and obtain the DGLAP evolution equation for the gluon structure 
function. If k ^ 1, we can use the asymptotic formula for Ei and obtain: 

xG.(.,Q^) = ^-^fj^j^^ '!^[c + ln(Ka{x',ri)} , 
where Ql{x') is the solution of the equation: 

In Fig.^ are plotted the contours of k for a nucleon target that give an idea in which kinematic 
region we expect big SC. 

3.2 Theory status of the Mueller formula. 

In this section we shall recall the main assumptions that have been made to obtain the 
Mueller formula. 

1. The gluon energy (x) should be high (small) enough to satisfy eq. (0) and as/n(l/x) < 
1. The last condition means that we have to assume the leading ln{l/x) approximation of 
perturbative QCD for the nucleon gluon structure function. 

2. The DGLAP evolution equations hold in the region of small x or, in other words, 
asln{l/r\) < 1. One of the lessons from HERA data is the fact that the CLAP evolution can 
describe the experimental data. 
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Figure 5: The contours of k for Nucleon, Ca and Au. 
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These two assumptions mean that we describe the gluon emission in so called Double Log 
Approximation ( DLA) of perturbative QCD, or in other words, we extract from each Feynman 
diagram of the order the contribution of the order (as Inl/a; InQ^/Qo neglecting all 
other contributions of the same diagram. In terms of the DGLAP evolution, we have to assume 
that the DGLAP evolution equations describe the gluon emission in the region of small x. 
However, the first assumption is very important for the whole picture, since it allows us to 
treat successive rescatterings as independent and simplifies all formulae reducing the problem 
to an eikonal picture of the classical propagation of a relativistic particle with high energy 
{E ^ where fi is the scattering radius in the nuclear matter) through the nucleus. The 
second one simplifies calculations but we can consider the BFKL evolution [|12] instead of the 
DGLAP one. 

3. Only the fastest partons {GG pairs) interact with the target. This assumption is an 
artifact of the Glauber approach, which looks strange in the parton picture of the interaction. 
Indeed, in the parton model we rather expect that all partons not only the fastest ones should 
interact with the target. In the next section we will show that corrections to the Glauber 
approach due to the interaction of slower partons are essential in QCD too. 

4. There are no correlations (interaction) between partons from the different parton cascades 
(see Fig. 4 ). This assumption means that even the interaction of the fastest GG-pair was taken 
into account in the Mueller formula only approximately and we have to assume that we are 
dealing with large number of colours to trust the Mueller formula. Indeed, it has been proven 
that correlations between partons from different parton cascades lead to corrections to the 
Mueller formula of the order of where Nc is the number of colours ( see Ref. and 
references therein for detail discussions on this subject). 

5. There are no correlations between different nucleons in a nucleus. 

6. The average for GG pair-nucleon interaction is much smaller than Ra- 

The last two are usual assumptions to treat nucleus scattering. We have used the specific 
Gaussian parameterization for bt dependence. Also, one can easily generalize our formula in 



more general Wood-Saxon parameterization 

3.4 The modified Mueller formula. 

The next step of our approach is to give an estimate of the SC using the Mueller formula. 
However, before doing so, we have to study how well works the DLA of perturbative QCD 
which was heavily used in the derivation of the Mueller formula. Let us recall that the solution 
of the DGLAP evolution equations can be easily found in the moments space. For any function 
f{x) we define the moment f{uj) as 

fiuj) = C dxx-fix) . 
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Note that the moment variable uj is chosen such that the = moment measures the number 
of partons, and the moment uj = \ measures their momentum. An alternative moment vari- 
able = — 1 is often found in the literature. The x-distribution can be reconstructed by 
considering the inverse Mellin transform, which for the gluon distribution reads: 

2 m Jc 

where the contour of integration C is taken to the right of all singularities and function gin is 
defined by the initial gluon distribution at = Qq. The anomalous dimension 'y{uj) has to be 
calculated in perturbative QCD and can be written in the form: 



asN, 1 , 2alN^a3) 1 
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7(^) = ^ ■ - + ■ - + O(^) + 0(as) . (20) 

IT U TT* CJ* UJ'^ 



In the DLA we take only the first term of this series, namely. 

In the BFKL evolution equation all terms of the order (^)" have to be taken into account. 
They generate the BFKL anomalous dimension of the form: 

where 'j^^^^{uj = ujl) = 1/2. The main qualitative property of the BFKL anomalous dimen- 
sion is the fact that it cannot exceed the value 1/2. 

The momentum conservation means that ''j{uj = 1) = 0. None of the DLA or the BFKL 
anomalous dimension satisfies this equation, because they give the good approximation to the 
full anomalous dimension only in the region of small values of uj or, in other words, in the region 
of small X. 




The DLA anomalous dimension leads to the simple evolution equation: 

■xG(x,g2). (21) 



d'^xG{x,Q'^ _ Oi^Nc 2\ 



91n(l/a;)91nQ2 7^ 

Now let us estimate how well works the DLA. In all our numerical estimates we use the GRV pa- 
rameterization 1^ for the nucleon gluon distribution, which describes all available experimental 
data quite well, including recent HERA data at low x. Moreover, GRV is suited for our purpose 
because (i) the initial virtuality for the GLAP evolution is small (Qg ~ 0.256*6^^) and we can 
discuss the contribution of the large distances in MF having some support from experimental 
data; (ii) in this parameterization the most essential contribution comes from the region where 
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aglnQ"^ ~ 1 and a^/nl/x ~ 1. This allows the use of the double leading log approximation of 
pQCD, where the MF is proven [^. It should be also stressed here, that we look at the GRV 
parameterization as a solution of the DGLAP evolution equations, disregarding how much of 
the SC has been taken into account in this parameterization in the form of the initial gluon 
distribution. 

However, in spite of the fact that the GLAP evolution in the GRV parameterization starts 
from very low virtuality ( Ql ~ 0.25GeV^) it turns out that the DLA still does not work quite 
well in the accessible kinematic region (Q^ > lGeV'^,x > 10^^). To illustrate this statement 
we plot in Fig. 5 the ratio: 



< 



> 



ain(l/a;)91n( 



TT 



This ratio is equal to 1 if the DLA holds. From Fig. 5 one can see that this ratio is rather 
around 1/2 even at large values of Q^. 



1.5 



1.0 



0.5 



0.0 



Q = 1 . Gev 

= 2.5 Gev' 

=5. Gev' 

=10. Gev' 

=20. Gev' 



1/ 



0.0 5.0 10.0 15.0 20.0 

y=ln(1/x) 



Figure 6: The ration for different values of in the GRV parameterization 



We can understand why the corrections to the DLA is so big modeling the complicated 
expression for 7 of eq. (^0|) by simple formula |]19| [|: 



-rM = ^-{^-l}. (22) 

TT U 

Eq. ( P^ ) has correct the DLA limit at small uj and it satisfies the momentum conservation 
7(0; = 1) = 0). The typical values of uj in all available parametrizations,even in the GRV 



^ We are very grateful to Yu. Dokshitzer for enlighting discussions on this problem during the RHIC'96 
Workshop 
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, which is the closest to the DLA, is < c<j > 0.5. Therefore, we have about 50% correction 
to the DLA. Therefore, the DLA cannot provide a rehable estimates for the gluon structure 
function. 

On the other hand, our master equation (see eq. (|T2p) is proven in DLA. Wilhng to develop 
a realistic approach in the region of not ultra small x {x > 10~^) we have to change our master 
equation ( eq. (|T^ ). We suggest to substitute the full DGLAP kernel ( the full expression of 
eq. (|20|) ) in the first term of the r.h.s. This procedure gives 



- I ll ^^x'G--(x',Q'^). (23) 



xGa{x,Q^) = xG'A(x,Q2)(eg. (12)) + AxG^«^(x,Q2 
asNc dx^dQ'^ 

TT Jx JQl X' Q 

The above equation includes also AxG'j^^ [x, Qq) as the initial condition for the gluon distribu- 
tion and gives AxG'^^ [x^ Q"^) as the first term of the expansion with respect to kg- Therefore, 
this equation is an attempt to include the full expression for the anomalous dimension for the 
scattering off each nucleon, while we use the DLA to take into account all SC. Our hope, which 
we will confirm by numerical calculation, is that the SC are small enough for x > 10~^ and 
we can be not so careful in the accuracy of their calculation in this kinematic region. Going to 
smaller x, the DLA becomes better and eq. (EBI) tends to our master equation 



The gluon structure function for nucleon (A = 1 ). 



In this subsection we are going to check how eq. (|23D describes the gluon structure function 



for a nucleon, which is our main ingredient in the Mueller formula. We calculate first the ratio 

^ _ xG^{x,Q^){eq. (23)) 



for A = 1, which is shown in Fig.[^. From this ratio we can see the general behavior of the SC 
as a function of ln{l/x) and Q^. In the region of the HERA data, 3 < ln{l/x) < 10, and 
> 2 GeV^^, the SC are not bigger than 15 %. The SC give a contribution bigger than 20 % 
only at very small value of x, where we have no experimental data. 

In the semiclassical approach (see 0), the nucleon structure function is supposed to have 
and x dependence as 

xG'^(x,Q2) oc {g2}<7> |1}<-> . (25) 

X 

We can calculate both exponents using the definitions 

gln(xG^(x,Q^)) 

< "-^ > = TTT-, — • (2d) 

(91n(l/x) ^ ^ 
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Figure 7: The SC for nucleon (A=l) as a function ofln{l/x) and , where ratio Ri compares 
xG^ with xG (GRV) distribution. 



< 7 > 



d\n{xG^{x,Q^)) 



(27) 



The eq.(P^) gives the average value of the effective power < u > of the gluon distribution, 
xG{x, Q^) oc which is suitable to study the small x behavior of the gluon distributions. 

Fig]^ shows the calculation of < a; > the nucleon distribution for eq. (^) and for GRV gluon 
distribution, both as functions of ln{l/x) for different values of Q^. From the figure, we can 
see that the effective powers of xG^{A = 1) and xG{GRV) have the same general behavior in 
the small x limit but the nucleon distribution is slightly suppressed. We calculate also, in the 
same kinematical region, the exponent < 7 >, given by eq (|27|) . This is the average value of 
the anomalous dimension, which describes the effective dependence of the distribution in 
variable. Figs.^ shows < 7 > for the nucleon and GRV distributions, indicating that the 
dependence is slightly soften by the SC. 

Comparing figures |^ and |^, we can conclude that even these more detailed characteristic of 
the gluon structure function have not been seriously affected by the SC in the nucleon case. 

We also use the DGLAP evolution equations to predict the value of the deep inelastic struc- 
ture function F2 from the xG"^ gluon distribution. Summing the DGLAP evolution equations 
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Figure 8: The effective power < cu > calculated for xG^{A = 1) and the GRV distribution. 
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Figure 9: The effective power of dependence calculated for xG^{A — 1) and the GRV 
distribution. 



16 



for each quark flavor, the function F2 may be written 




where the sea quark distributions have been neglected in comparison with the gluon distribution. 
Figjl^ shows the prediction for F2 from xG^ and from the GRV distribution, compared with 
experimental data. As we can see, the magnitude of the suppression due to the SC is less than 
10% in the region of the HERA data and this suppression is smaller than the experimental 
error. 




Figure 10: F2 from xG^ and the GRV distribution, compared with experimental data J^. 

From the above results we can conclude that eq. (^) gives a good description for the gluon 
structure function for nucleon and describes the available experimental data. Therefore, it can 
be taken as a correct first approximation in the approach to the nucleus case. 

3.6 The gluon structure function for nucleus. 

In the framework of perturbative approach it is only possible to calculate the behavior of 
the gluon distribution at small distances. The initial gluon distribution should be taken from 
the experiment. Actually the initial virtuality Ql should be big enough to guarantee that we 
are dealing with the leading twist contribution. Our main assumption is that we start the 
QCD evolution with a small value of Qq considering that the MF is a good model for high twist 
contributions in DIS off nucleus. 
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Table 1: Values of Rin and a for parameterization Ri = Rin A 







Q' = 


IGeV^ 


Q' = 


lOGeV^ 


X 




RlN 


a 


RlN 


a 


10- 


-■2 


0.94 


0.0416 


0.98 


0.014 


10- 


-3 


0.92 


0.0616 


0.94 


0.034 


10- 


-4 


0.88 


0.094 


0.92 


0.0563 


10- 


-5 


0.8 


0.145 


0.86 


0.093 



The scale of the SC governs by the value of k^, namely they are big for ka > I and small for 
ka < I. Fig.4 shows the plot of = 1 for different nuclei. One can see that the SC should be 
essential for heavy nuclei starting from Ca at the accessible experimentally kinematic region. 

Now we extend the definition of Ri for the nucleus case 

AxG%^''{x,Q^) ' ^^""^ 

where the numerator is calculated using eq.(^). Figure [TTl shows the results for the calculations 
of Ri as a function of the variables ln{l/x), InQ"^ and A^/^. Fig.^a presents the ratio Ri for 



two different values of and for different nuclei. The suppression due to the SC increases with 
ln{l/x) and is much bigger than for the nucleon case. For A = 40 (Ca) and = lOGeV"^, 
the suppression varies from 4% for ln{l/x) = 3 to 25 % for ln{l/x) = 10. For A = 197 (Au) 
the suppression is still bigger, going from 6% to 35% in the same kinematic region. Fig. |TT|b 
shows the same ratio for different values of for the gold. The suppression decreases with 
Q^. Figs. |TT]c and d show the Ri ratio as a function of A^^^ and x for a fixed value of Q^. 
As expected, the SC increases with A. An interesting feature of this figure is the fact that the 
curves tend to straight lines as x increases. It occurs because, as x grows, the structure function 
xG{GRV) becomes smaller, and the correction term of (^) proportional to k, dominates. Since 
K is proportional to A^^^, the curves behave as straight lines. The decrease of suppression with 
is illustrated in more detail in Figs. |Tl|e and f which presents Ri as a function of In for 
different values of x for Ca and Au, respectively. The effect is pronounced for small and x 
and diminishes as InQ^ increases. 

This picture ( Fig.lTl] ) shows also that the gluon structure function is far away from the 



asymptotic one. The asymptotic behavior i?i — > 1 ( see Figs. lie and f ) occurs only at very 



high value of as well as in the GLR approach ( see ref. [ESI ). The asymptotic A-dependence 



( Ri oc A"^) ) has not been seen in the accessible kinematic range of and x ( see Figs. [Tl| c 
and d and Table 1 ). This result also has been predicted in the GLR approach We want also 
to mention that parameterization Ri = Rin A~°' does not fit the result of calculations quite 
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Figure 11: Ri as a function of ln{\/x), InQ^ and A^^^: a) Ri as a function of ln{\/x) for 
different nucleus and different values of A; b) Ri as a function of InQ'^ for different values of 
xb for Au; c) and d) Ri as a function of A^/^ for different Q^; e) and f) Ri dependence on 
for Ca and Au. 
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well for IGev'^ < < 20GeV'^ and 10"^ < x < 10^^. For x ~ 10"^ the parameterization 
Ri = i?i7v — R' As with parameters Rin and R' for each value of Q^, works much better 
reflecting that only the first correction to the Born term is essential in the Mueller formula. 

We extend also the calculation of the exponents < u > and < 7 > of the semiclassical 
approach for the nuclear case. We calculate the effective power of the nuclear gluon distribution 
< u > using the expression 

dHi/x) ■ 

Fig]l^ shows the results as functions of ln{l/x) for different values of and different nucleus. 
The SC decreases the effective power of the nuclear distribution, giving rise to a flattening of 
the distribution in the small x region. 
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Figure 12: < a; > for different values of and A. 

It is also interesting to notice that at small values of Q^, the effective power tends to be 
rather small, even in the nucleon case, at very small x. However it should be stressed that 



the effective power remains bigger than the intercept of the so called "soft" Pomeron pi 
even in the case of a sufficiently heavy nucleus (Au), for > IGeV"^. Nowadays, many 
parameterizations p2[ with matching of "soft" and "hard" Pomeron have appeared triggered 
by new HERA data on diffraction dissociation . These parameterization used Pomeron-like 
behavior namely, xG{x,Q'^) oc x~^^^ \ However, if the Pomeron is a Regge pole, u cannot 
depend on Q"^, and the only reasonable explanation is to describe uj^Q"^) as the result of the SC. 
Looking at Fig.|T^ we can claim the SC from the MF cannot provide sufficiently strong SC to 
reduce the value of u to 0.08, a typical value for the soft Pomeron [^], at least for Q'^ > IGeV"^. 

The calculation of the effective value of the anomalous dimension 7 may help us to estimate 
what distances work in the SC corrections. This effective exponent is given by 



< 7 > 



(91n(xG'^(x,g2)) 



(31) 
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Figure 13: < 7 > for different Q"^ and A. 
Fig. O shows the resuhs as functions of ln{l/x) for different values of and for two nuclei. 



We see that the values of 7 at ln{l/x) < 5, for both Ca and Au, is very close to the results 
for GRV and for nucleon case. At smaller values of x, the anomalous dimension presents a 
sizeable reduction, which increases with A. For ln{l/x) > 15, < 7 > tends to zero unlike in 
the DGLAP evolution equations ( see Fig. 10 for the GRV parameterization). Analysing the 
dependence, we see that < 7 > is bigger than 1 only for = 0.5 GeV^. For = 1.0 GeV^, 
the anomalous dimension is close to 1/2, and for > 5.0 GeV^ it is always smaller than 1/2. 

Using semiclassical approach, we see that 

'^oc^(Q2)\ (32) 

and if 7 > 1, the integral over rj in the master equation (^) becomes divergent, concentrating 
at small distances. 

If 1 > 7 > 1/2, only the first SC term, namely, the second term in expansion of the master 
equation, is concentrated at small distances, while higher order SC are still sensitive to small r± 
behavior. Fig.O shows that this situation occurs for > 1 GeV^, and even for = IGeV^ 



at very small values of x. We will return to discussion of these properties of the anomalous 
dimension behavior in the next section. 

3. 7 The gluon life time cutoff. 

In the DIS the incident electron penetrates the nucleus and radiates the virtual photon 
whose lifetime t^* oc — 12711. We can recover three different kinematic regions: 



^- '^■y* ~ mx ^ Rnn, where Rnn is the characteristic distances between the nucleons 
of the nucleus. This virtual photon can be absorbed only by one nucleon and the total cross 
section is cr(7*A) = Aa{Yp)- 

2. Ra > '^7* = ^ > Rnn, where Ra is the nucleus radius. In this kinematic region the 
virtual photon can interact with the group of nucleons. However, a{'j*A) is still proportional 
to A since the number of nucleons in a group is much less than A. 
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3. T^* = — > Ra- Here, before reaching the front surface of the nucleus, the virtual 

' mx ^ I o 7 

photon "decomposes" in the developed parton cascade which then interacts with the nucleus. 
It can be shown PB[ that the absorption cross section of the virtual photon will now be propor- 

2 

tional to the surface area of the nucleus cr(7*y4) oc As, because the wee partons of the parton 
cascade are absorbed at the surface and do not penetrate into the centre of the nucleus. 

Everything that we have discussed have been calculated in the third kinematic region. 
For the RHIC energies we have to develop some technique how to penetrate into the second 
one. To do this we have to remember that the opacity ( or ) actually depends on the 
longitudinal part of the momentum transfer ( ) which could be calculated in terms of x and 
x' of our master equation (|1^) , namely, = (x + x')m ( see Ref.[jl[ ). Recalling that opacity 
VL oc r\xG{x,qz,-^) S(bt,qz) we see that q^ - dependence enters to two factors: to gluon 
structure function and to the nucleon profile function. We know how to take into account the 
qz - dependence of the gluon structure function ( see Ref. where the DGLAP equation for 

7^ is written ). However we neglected this effect in our present estimates, hoping that this 
dependence occurs on the hadron scale and cannot change too much the dependence of the SC 
on the number of collisions during the passage through the nucleus. 

The dependence of the profile function S{bt, qz) on qz have been discussed and in the Gaus- 
sian parametrization it can be factor out in the form: 



2 

A „2 



Sibt,qz) = Sibt) ■ L{qz) = Sibt) ■ e--^^ . (33) 

This qz - dependence takes into account the fact that the virtual gluon can interact with the 
target only during the finite time r = 1/mx undergoing pr < pRa collisions. Using eq. (p3D, 
we can obtain: 

xGa{x,Q') = AxGr,{x,Q') - (' [ L{qz) x' Gm{x' ,Q") (34) 

TT Jx J X' U'^ 



^ dx^ r^^ drj ln(L(g.)«:c(x',r,^))+ E,{L{qz) ko{x' , r^) )} 

ir^ Jx x' J4^ rj 



Fig. |T4| shows the result of our calculations. Comparing Fig.|TT| with this picture, one can see 
that the finite life time of the virtual gluon affects the behavior of the gluon structure function 
only at sufficiently large x { x > 10~^ ) diminishing the value of the SC in this kinematic 
region. This effect turns out to be very important for the RHIC energies and has to be studied 
in more details. 



4 First corrections to the Glauber ( Mueller) Approach. 

In this section we discuss the corrections to the Glauber approach (the Mueller formula of 
eq. ( p!2D . To understand how big could be the corrections to the Glauber approach we calculate 
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the second iteration of the Mueller formula of eq. (0). As has been discussed,eq. describes 
the rescattering of the fastest gluon ( gluon - gluon pair ) during the passage through a nucleus 
( see Fig.l ). In the second iteration we take into account also the rescattering of the next to 
the fastest gluon. This is a well defined task due to the strong ordering in the parton fractions 
of energy in the parton cascade in leading ln{l/x) approximation of pQCD that we are dealing 
with. Namely: 

xb < Xn < ■■■ < xi < 1 ; (35) 
where 1 corresponds to the fastest parton in the cascade. 

Therefore, in the second interaction we include the rescatterings of the gluons with the 
energy fraction 1 and Xi ( see Fig.[l^ ). Doing the first iteration we insert in eq. (|l^) Gn{x, Q^) = 
G%^^ {x, Q^). For the second iteration we calculate the gluon structure function using eq. ([l2[) 
substituting 

xGj, = ^^M^ _ xG^^^(x,Q^) ; (36) 



where xG\ is the result of the first iteration of eq. (|T2|) that has been discussed in details in 
section 3. 



Fig.|T6| shows the need to subtract xG%^^ in eq. ( p6l) making the second iteration. Indeed, 



in the second iteration we take into account the rescattering of gluon 1' - gluon 2' pair off 



a nucleus. We picture in Fig.[l^ the first term of such iteration in which G11G2' pair has no 
rescatterings. It is obvious that it has been taken into account in our first iteration, so we have 
to subtract it to avoid a double counting. 
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Figure 15: The interaction with nucleons that is taken into account in the second iteration 
M ueller formula. 
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Figure 17: Second iteration calculations for Ri, < uj >, and < 7 > for Ca and Au. 
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One can see in Figs-O that the second iteration gives a big effect and changes crucially 



< 7 >, and < uo >. The most remarkable feature is the crucial change of the value of the 
effective power uj{Q'^) for the "Pomeron" intercept which tends to zero at HERA kinematic 
region, making possible the matching with "soft" high energy phenomenology. It is also very 
instructive to see how the second iteration makes more pronounced all properties of the behavior 
of the anomalous dimension ( < 7 >) that we have discussed. The main conclusions which we 
can make from Figs.|T^ are: (i) the second iteration gives a sizable contribution in the region 
X < 10~^ and for x < 10^^ it becomes of the order of the first iteration; (ii) for x < 10~^ 
we have to calculate the next iteration. It means that for such small x we have to develop a 
different technique to take into account rescatterings of all the partons in the parton cascade 
which will be more efficient than the simple iteration procedure for eq. (|12]). However, let us 
first understand why the second iteration becomes essential to establish small parameters that 
enter to our problem. 

As has been discussed, we use the GLAP evolution equations for gluon structure function 
in the region of small x. It means, that we sum the Feynman diagrams in pQCD using the 
following set of parameters: 



1 Q'^ I 

as < 1 ; as In - < 1 ; as In —2 < 1 ; as In —2 In - ^ 1 . (37) 

X C^Q C^Q X 



The idea of the theoretical approach of rescattering that has been formulated in the GLR paper 
[H is to introduce a new parameter []: 



xG{x,Q') (38) 



A 



and sum all Feynman diagrams using the set of eq. ( P7| ) and k as parameters of the problem, 
neglecting all contributions of the order of: as, as as ln(l/x), as ln(l/a;) k, as lia{Q'^/Ql) and 
as hi{Q'^/Ql) K. It should be stressed that Mueller formula gives a solution for such approach. 
Indeed, eq. (^2]) depends only on k absorbing all 

(as ln(Q^/(5o) ln(l/a;))'^ contributions in xG{x,Q'^). However, it is not a complete solution. 
To illustrate this point let us compare the value of the second term of the expansion of eq. (|T2|) 
with respect to (j{rf) with the first correction due to the second iteration in the first term of 
such an expansion. In other words we wish to compare the values of the diagrams in Fig.lTB| b 



and Fig.|T8|a. The contribution of the diagram of Fig.p!8|a is equal: 

AxG{x,Q'){Fig.l%a) = / ^ / dQ'^ n\x\'^) , (39) 

2tt^ J x' J 4 

where x' and Q'^ are the fraction of energy and the virtuality of gluon 1 in Fig.[T^a 

^ In the GLR paper the notation for k was W, but in this paper we use k to avoid a misunderstanding since, 
in DIS, W is the energy of interaction. 
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The diagram of Fig.lT8|b contains one more gluon and its contribution is: 



K [X , 



TT^2 n J x' Q'2 J 
oc ln(l/x) IniQ^Ql) AxG{x,Q^)iFig.l8a)) , (40) 



where x' {x") and Q'^ ( Q"^) are the fraction of energy and the virtuahty of gluon 1 (1') 
respectively in Fig.|18|b. Therefore, eq. (^) gives the contribution which is of the order of 
eq. ( P^D in the kinematic region where the set of parameters of eq. ( P7D holds. It means also 
that we need to sum all diagrams of Fig.p!8|b type to obtain the full answer. In the diagram of 
FigjT^ not only one but many gluons can be emitted. Such emission leads to so called "triple 
ladder" interaction, pictured in Fig.|18|c ( see ref.0| ). This diagram is the first from so called 
"fan" diagrams of Fig. |l8]d. To sum them all we can neglect the third term in eq. (0) and 
treat the remained terms as an equation for xGlXjQ"^). It is easy to recognize that we obtain 
the GLR equation 0[|lOl- Generally speaking the GLR equation sums the most important 
diagrams in the kinematic region where as ln(l/x) ln{Q'^/Ql) ^ 1 and k < 1. 



5 The general approach. 

5.1 Why equation? 

We would like to suggest a new approach based on the new evolution equation to sum 
all SC. However, first of all we want to argue why an equation is better than any iteration 
procedure. To illustrate this point of view let us differentiate the Mueller formula with respect 
to y = ln(l/x) and C, = InQ^. It is easy to see that this derivative is equal to 

d^xGix,Q^) _ 4 _ ^4,) 



dy TT^ 

The nice property of eq. (^) is that everything enters at small distances, therefore everything 
is under theoretical control. Of course, we cannot get rid of our problems changing the pro- 
cedure of solution. Indeed, the nonperturbative effects coming from the large distances are 
still important but they are all hidden in the boundary and initial conditions to the equation. 
Therefore, an equation is a good ( correct ) way to separate what we know ( small distance 
contribution) from what we don't ( large distance contribution). 

5.2 The generalized evolution equation. 

We suggest the following way to take into account the interaction of all partons in a parton 
cascade with the target. Let us differentiate the Mueller formula over y = ln(l/x) and ^ = 
HQ^Ql). It gives: 

d^xGA{y,0 ^ 2R\Q^ 
dydi vr 



2 {G + hiK + Ei{k) } . (42) 
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Figure 18: Corrections to the Glauber approach. 
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Rewriting eq. (^) in terms of n given by 

^ = ^^/Ga{x,Q') (43) 



we obtain: 

+ — 7^ = I <^ + In/t + Ei{K{y,0) } = FW ■ 44 

Now, let us consider the expression of eq. as the equation for k This equation has the 
following nice properties: 

1. It sums all contributions of the order [asV^)"' absorbing them in xGA{y,C,), as well as 
all contributions of the order of k". Therefore, this equation solves the old problem, formulated 
in Ref.0] and for ^ oo eq. gives the complete solution to our problem, summing all 
SC; 

2 .The solution of this equation matches with the solution of the DGLAP evolution equation 
in the DLA of perturbative QCD at k ^ 0; 

3. At small values of /t ( k < 1 ) eq. (|^) gives the GLR equation. Indeed, for small k we 
can expand the r.h.s of eq. (|^) keeping only the second term. Rewriting the equation through 
the gluon structure function we have 



2 



xG{x,Q') - ^^{xG{x,Q')Y , (45) 



which is the GLR equation with the coefficient in front of the second term calculated by 
Mueller and Qiu 0. 



4. For a^yS, ~ 1 this equation gives the Glauber ( Mueller ) formula, that we have discussed 
in details. 



5. This equation almost coincide with the equation that L.Mclerran with collaborators |29 
derived from quite different approach and with different technique. We are sure that almost 
will disappear when they will do more careful averaging over transverse distances. 

Therefore, the great advantage of this equation in comparison with the GLR one is the fact 
that it describes the region of large k and provides the correct matching both with the GLR 
equation and with the Glauber ( Mueller ) formula. 



Eq. ( |41D is the second order differential equation in partial derivatives and we need two 
initial ( boundary ) conditions to specify the solution. The first one is obvious, namely, at fixed 

y and Q"^ ^ oo 



2Q''R\ 
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The second one we can fix in the following way: aX x = xq {y = yo) which is small, namely, in 
the kinematic region where asy^ < 1 



= 2Q^]i^ xGAix,Q^) , (46) 



where xGa is given by the Mueller formula ( see eq. ([121)). Practically, we can take Xq = 10 ^, 
because corrections to the MF are small at this value of x = Xq. 

5.3 The asymptotic solution. 

First observation is the fact that eq. ( PI) has a solution which depends only on y. Indeed, 
one can check that n = i^asympiy) is the solution of the following equation: 



-^( ^asymp ) • l^"^) 



The solution to the above equation is: 

rKasymp(y) dl^' 



y-yo. (48) 



It is easy to find the behavior of the solution to eq. (^81) at large value of y since F(k) — »• 
as In K' at large k { as = ^ ces )■ It gives 

Kasymp ^ "sZ/ ln(«s?/) at ag?/ > 1 • (49) 

At small value of y, F{k) asi^ and we have: 



asymp ' "asymp 



f^asympiV = Z/o) C^^^^"^") . (50) 



The solution is given in Fig.|T9| for as = 1/4 in the whole region of y for different nuclei in 
comparison with our calculations based on the MF. We chose the value of i^asympiy = 2/o) from 
eq. (|46D . We claim this solution is the asymptotic solution to eq. (0) and will argue on this 
point a bit later. 

For nuclei the SC incorporated in the asymptotic solution turn out to be much stronger 
than the SC in the Glauber approach for any > IGeV'^ at a; > 10~^. In this kinematic 
region the solution of eq. (H3) is drastically different from the Glauber one. 



A general conclusion for Fig.lTUI is very simple: the amount of shadowing which was taken 



into account in the MF is not enough , at least for the gluon structure function in nuclei at 
X < 10~^ and we have to solve eq. (^4|) to obtain the correct behavior of the gluon structure 
function for nuclei. 
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Figure 19: The Glauber approach and asymptotic solution for different nuclei. 

Now, we would like to show that solution eq. (|^) is the asymptotic solution of the new 
evolution equation. In order to check this we need to prove that this solution is stable. It 
means that if add a small function Ak and searching for the solution to the equation in the 
form K = Kasymp + ^1^, we have to prove that Ak turns out to be small. namely, A/t ^ k. 
The following linear equation can be written for Ak: 

d'AKiy,^ dAKjy,^ _ dF{K) 

dydi ^ dy ~ dK 1-=-— (^) ^'^^^'^^ ^^^^ 

In Ref.[|l| was proven, that the solution of eq. (pTD is much smaller than k. 



Therefore the asymptotic solution has a chance to be the solution of our equation in the 
region of very small x. To prove that the asymptotic solution is the solution to the equation 
we need to solve our equation in the wide kinematic region starting with our initial condition. 
We managed to do this only in semiclassical approach. 

5.4 Semiclassical Approach. 

The semiclassical approach has been adjusted to the solution of the nonlinear equation of 
eq.(0)-type in Refs. 0, [S^, |^ ( for simplicity, we assume that as is fixed ). 

In the semiclassical approach we are looking for the solution of eq. (0) in the form 

K = e^ (52) 

where S" is a function with partial derivatives: || = a; and f| = 7 which are smooth function 
of y and ^. It means that 

d^S dS dS 

< ■ TTT = (53) 



d^dy dy 
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Using eg. (^31) , one can easily rewrite eq.(^) in the form 

dS dS dS 
dy dy 



^^^^ + ^ = e-'Fie')^^S) (54) 



or 



cu(7 + l) = $(5) (55) 

We are going to use the method of characteristics ( see, for example, ref.|3^). For equation 
in the form 

F{^,y,S,^,uj) = (56) 

we can introduce the set of characteristic lines {^{y), S{y),u!{y),'y{y)), which satisfy a set of 
well defined equations (see, for example, Refs. |^ |^ for the method and Ref 0] for de- 
tailed calculation). Using eq. (^4]) and eq. (|55|) , we obtain the following set of equations for the 
characteristics: 

dy (7+1)=' df, (7 + 1)'^'' dy ^7 + 1' ^ ' 

where = The initial condition for this set of equations we derive from eq.(p6D, namely 



7o 



So = lnKin{yo,^o 
dlnKin{yo,0 



(58) 
Co 



The main properties of these equations have been considered in Ref. [0] analytically, however, 
here, we restrict ourselves mostly the numeric solution of these equations. 

We set the initial condition y = yo = 4.6 [xb = 10~^), where the shadowing correction is 
not big and the evolution starts from 7 < 0. In this case d'~f/dy>0 and the value of 7 increases 
. At the same time dS/dy < and 5* decreases if 70 < —1/2. With the decrease of S, the value 
of becomes smaller and after short evolution the trajectories of the nonlinear equation start 
to approach the trajectories of the DGLAP equations. We face this situation for any trajectory 
with 7o close to -1. If the value of 70 is smaller than — ^ but the value of 5*0 is sufficiently 
big, the decrease of 5* due to evolution cannot provide a small value for $'(5*) and 7 increases 
until its value becomes bigger than — ^ at some value of y = yc- In this case for y > yc the 
trajectories behave as in the case with 70 > — |. For 70 > —1/2, the picture changes crucially. 
In this case, dS/dy > , d'-f/dy > and both increase. Such trajectories go apart from the 
trajectories of the DGLAP equation and nonlinear effects play more and more important role 
with increasing y. These trajectories approach the asymptotic solution very quickly. 
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For the numerical solution we use the 4th order Runge - Kutta method to solve our set of 



equations with the initial distributions of eq. (|58D . The result of the solution is given in Figs.^ 
and pri In these figures we plot the bunch of the trajectories with different initial conditions. For 



the nucleon ( Fig.^ ) we show also the dependence of 7 along these trajectories. One can notice 
that the trajectories behave in the way which we have discussed in our qualitative analysis. It is 
interesting to notice that the trajectories, which are different from the trajectories of the GLAP 
evolution equations, start at y = = 4.6 with the values of between O.bGeV^ and 2.56*6^^ 
for a nucleon. It means that, guessing which is the boundary condition at = Qq = 2.5GeV^, 
we can hope that the linear evolution equations ( the DGLAP equations) will describe the 
evolution of the deep inelastic structure function in the limited but sufficiently wide range of 

In Figs. ^ and|21| we plot also the lines with definite value of the ratio R = ^ '^(^ ' '^^^^^'^a "'g^ j^p)""*^""^ 
(horizontal lines). These lines give the way to estimate how big are the SC. One can see that 
they are rather big. 

Contour plot for Nonlinear eq. 
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Figure 20: The trajectories and contour plot for the solution of the generalized evolution 

PmmHnn fnr N R — a:G(3;,Q^)(generafi2ed equation) 

equauon joi i\. n - xG{x,q^){glap) 

We have discussed only the solution with fixed coupling constant which we put equal to 
as = 0.25 in the numerical calculation. The problem how to solve the equation with running 
coupling constant is still open. 

5.5 The generalized evolution equation versus the GLR equation. 

In Ref. [0 we studied in detail the solution to the GLR equation in the same semiclassical 
approximation. Our conclusion is that the GLR equation gives much stronger SC than the 
generalized evolution equation. This difference we can see comparing the solution to the both 
equation in the region ultra small x. 
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Figure 21: The trajectories and contour plot for the solution of the generalized evolution equation 
for Ca and Au. R = "^^"''^I^g'ja)"^;)""*'""^ 

Indeed, our asymptotic solution turns out to be quite different from the GLR one. Tlie 
GLR solution in the region of very small x leads to saturation of the gluon density ^ . 

Saturation means that k tends to a constant in the region of small x. The solutions of eq. (^if ) 
approach the asymptotic solution at a; — > 0, which does not depend on Q^, but exhibits 
sufficiently strong dependence of k on x ( see Fig.|l9| ), namely n oc as ln(l/x) lnln(l/x). The 
absence of saturation does not contradict any physics since gluons are bosons and it is possible 
to have a lot of bosons in the same cell of the phase space. We should admit that A. Mueller 
first came to the same conclusion using his formula in Ref.||^. 



6 Next steps. 

Here, we list our problems that have to be solved to complete our study of the SC : 

1. Calculation of F^{x, Q"^) to compare our calculation of the SC with the available exper- 
imental data. 

2. Recalculation of the SC using more reliable Wood-Saxon parameterization for profile 
function S{ht) instead of the Gaussian one. The form of the profile function especially essential 
to obtain a reliable estimates for the SC in the region of the moderate x < 2mR^ ■ 

3. Solution of the generalized evolution equation for running as- The experience of solving 
the GLR equation tells us that there is a principal difference in the solutions for fixed and 
running as, namely, the critical line of the GRL equation appears only for running as 0. 
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We think that it is very important to study the generahzed equation with running as and to 
compare this solution with the solution of the GRL equation. 

4. We have discussed that for RHIC energies it is very important to study in more details 
the effect of the final life-time of th gluon in a nucleus. We plan to recalculate the SC replacing 
xG{x, Q"^) in our formulae by xG{x, Q^, Qz) for which the kernels of the evolution equations 
have been calculated in Ref. p|. 

5. In all our calculations we neglected the parton interaction inside GG + N scattering. Our 
estimates, which have been presented in section 2, shows that this interaction should be very 
important. Indeed, for example, in the Mueller formula we have to change the parameter kq 
due to the parton interaction inside the nucleon. This change is simple, the only that we need 
to do is to replace the number of collisions A/7cR\ by 



A 



+ 



7iR\ 7iR\ nR% 

in the definition of hg in eq. (p!2D . It means that all results will be the same but nucleus with 
the new effective number of nucleons: 

^eff - ^ + d2 ' 

where Ra = RqA^. Using our estimates for i?^ = ^GeV'"^ we can see that effective A for 
1 1 

the gold is Al^j = 9.6 instead Aa =6. For light nuclei the change is even more essential. 
Therefore, we are planning to take into account the parton interaction inside a nucleon as soon 
as possible. 

6. We have neglected all correlations between partons of the order which could be sizable 
in the case of the nucleus DIS. We suppose to study this problem using the technique that has 
been developed in Ref. p3 . 



7. Everywhere through the paper we used the DLA of perturbative QCD. However, the key 
assumption that simplify our theoretical approach was the aslB.{l/x) ~ 1 approximation. We 
plan to develop our approach in the case of the BFKL dynamic and, therefore, to get rid of 
our assumption that as Ib-IQ'^/Qq) ^ 1. We consider this generalization as an important step, 
since our result that we have no saturation of the gluon density in nuclei even at ultra small x 
could be an artifact our double log approximation of perturbative QCD. 



7 And what? 

We presented here our approach to the SC and a natural question arises:and what? What and 
how we can do for the RHIC physics. How our approach can help in creating of the reliable 
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Figure 22: Four stages of nucleus - nucleus collision. 

Monte Carlo code for nucleus - nucleus interaction at high energies. We are going to answer 
these hot question in this section. 

Let us consider first the space time structure of the nucleus - nucleus interaction ( see 
Figj2^). One can see four stages of this process: 

1. For time smaller that t^, where to is the time of the first parton - parton interaction, we 
have a very coherent system of parton, confined in our both nuclei. We know almost nothing 
about this system. 

2. At time to the first parton - parton interaction occurs and we believe that this interaction 
destroys the coherence of our parton system at the very instant. 

3. During time from to till t/j, where t/j is the hadronization time, we have a quark - gluon 
stage of the process. We believe that we can reach a simple and economic understanding this 
stage in framework of QCD. We also believe that new collective phenomena could be created 
in the nucleus - nucleus interaction during this stage of the process such as the Quark - Gluon 
Plasma mostly because of the high density of the produced gluons. For this stage we have the 
Monte Carlo codes based on QCD, the lattice calculation and a lot of beautiful ideas that has 
been discuss at this conference. 

4. The last stage - hadronization is a black box. Nothing is known, but the success of the 
Local- Hadron-Parton Duality in the description of the LEP data allows us to hope that this 
stage could be not very important for our understanding of the nucleus - nucleus collisions. 

Our approach can define the initial condition at t = to for the third stage. What can we 

do? 



1. We are able to calculate the inclusive cross section for gluons at t = to or, in other words, 
define the gluon distribution at t = to. Actually, it has been done by C.Escola and his 
collaborator and has been presented at this conference. We can only improve his treatment of 
the SC which was based on the GLR equation. However, let us discuss briefly the formula for 



the inclusive gluon cross section. It can be written using the factorization theorem ||36| in the 
form: 



ay api J pi 
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where the last factor is the hard gluon - gluon cross section and y and pt are rapidity and 
transverse momentum of produced gluon, respectively. One can see that this cross section is 
infrared unstable and diverges at small values of Pf The SC provides a natural scale that cut 
off this divergence. A rough estimate for this new scale can be done from equation 



1 



(see Fig. 4 ). For < Q){x) the gluon structure function xG{x,pf) oc pf and one can see that 
the number of gluon with transverse momenta smaller than pt = Qo{x) turns out to be very 
small. 

2. We can calculate also the double inclusive cross section which gives the two gluon 
correlation function at t = to- We would like to stress that for nucleus - nucleus collision this 
correlation function is big and have to be taken into account. Indeed, we have two different 



contribution to the double inclusive process, pictured in Fig.^: the production of two gluons 
from one parton cascade (see Fig.|23|a) and from two parton cascades ( see Fig.p3|b ) . However, 
for nucleus - nucleus collisions the first contribution is proportional to ^41^42 ( without the SC) 

while the second is much bigger and it is of the order of AfA^ ( without the SC and for 

the Gaussian profile function). Using our approach we can calculate the two gluon correlation 
function within better acuraccy than the above simple estimates. We hope, that these two 
observables: gluon distribution and two gluon correlation function will be enough for reliable 
description of the initial condition for the QCD motivated cascade during the third stage of 
our process. 

3. We think that these two observables: multiplicity of gluons and two gluon correlation 
will be enough to define the initial condition for current Monte Carlo codes. However, we think 



that these codes are doing something wrong. Indeed, we learned from A. Mueller [|T7l that 
correct degrees of freedom for parton cascading looks in the simplest way and which could be 
used for a probabilistic interpretation and therefore, they are natural degrees of freedom for 
Monte Carlo simulations are not quark and gluons but colourless quark - antiquark dipoles. 
The gluon structure function is the probabihty to find a colourless dipole with the size r± > 
Therefore, we think that the code should be written for such dipoles and their interaction. We 
shall answer the questions: (i) how to calculated the average multiplicity of dipoles with the size 
r_i_ and (ii) how to calculate the correlations between such dipoles. We are going to do this in 
the nearest future. 

4. Now we want to discuss a hot question how to mix the "soft" and "hard" Pomerons. 
The common way of doing such a mixture is to use the Glauber formula and replace in this 
formula (j{rj_ — > agoft + CThard- We think this is a correct procedure to obtain an estimate how 
important soft or/and hard processes. In section 3 we argued that this is the most economic 
way of doing which satisfies the s-channel unitarity. However, all Monte Carlo programs that 
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Figure 23: Double inclusive cross section in QCD. 
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(4) 



incx^S 



Figure 24: Total cross section in ag and orders of perturbative QCD. 



we know use for the calculation of (Jhard the factorization formula, namely 



^hard 



dpi xiG{x2,pI)x2G{x2, pI 



Pt 



which describes really the inclusive production of gluons. The factor 1/2 in front does not 
help because to find auard we need to calculate the real multiplicity but not the number of 
gluon line in the Feynman diagram. In Fig.|2^ we picture the corrections to the hard cross 
section considering the scattering of two mesons made from heavy quarks. Perturbative QCD 
is certainly a good tool to study such processes. From this picture one sees that including the 
inclusive cross section in the place of the total we missed the radiative correction to the the 
partial cross section with four quarks i the final state. 

Our way of doing is the following. We will write the Mueller formula or our more sofisticated 
approach for dipole (with size r± scattering with a nucleus. To find the proton - nucleus cross 
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section we need to calculate the integral: 



1 3 72 ik 

To find the wave function of the nucleon wc have to use a model, for example the constituent 
quark model or instanton liquid model. The nice feature of this formula that the typical 
will be of the order 1 GeV^^ due to the SC. It means that wc need to know the wave function 
at sufficiently small distances where we have some control from lattice calculations and QCD 
sum rules. This formula takes into account correctly hard process and give the factorization 
formula for the inclusive production. We suppose to do an estimates using the model for the 
nucleon wave function. If they will show that we need some admixture of the soft processes we 
will add to o-jy in the above formula in an usual phenomenologic way, using the model of, so 
called, soft Pomeron. 



8 Conclusions. 

We have two conclusions: 

1. We hope that we convinced you that we are on the way from our Really Highly Inefficient 
Calculation to your RHIC. Much work is need to clarify the initial condition for the QCD phase 
of nucleus - nucleus interaction and this is the first and the most important task which wc need 
to attack, since it will determine the correct degrees of freedom for further evolution of QCD 

cascades. 

2. Everything that we have talked about satisfies the third law of theoretical physics:^ny 
model is a theory which we apply to a kinematic region, where we cannot prove that this theory 
is wrong. We firmly believe that correct SC will provide the picture of the nucleus - nucleus 
interaction in which hard and semihard processes will play a crucial role with only small if any 
contamination of the soft contribution. 
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